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ABSTRACT 

Sky temperature map of the cosmic microwave background (CMB) is one of 
the premier probes of cosmology. To minimize instrumentally induced system- 
atic errors, CMB anisotropy experiments measure temperature differences across 
the sky using paires of horn antennas with a fixed separation angle, temperature 
maps are recovered from temperature differences obtained in sky survey through a 
map-making procedure. The instrument noise, inhomogeneities of the sky cover- 
age and sky temperature inevitably produce statistical and systematical errors in 
recovered temperature maps. We show in this paper that observation-dependent 
noise and systematic temperature distortion contained in released Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) CMB maps are remarkable. These errors 
can contribute to large-scale anomalies detected in WMAP maps and distort the 
angular power spectrum as well. It is needed to remake temperature maps from 
original WMAP differential data with modified map-making procedure to avoid 
observation-dependent noise and systematic distortion in recovered maps. 

Subject headings: cosmic microwave background — cosmology: observations - 
methods: statistical 



1. Map- making 

The COBE and WMAP missions measure temperature differences between sky points 



using differential rad iometers consisting of plus-horn and minus-horn (ISmoot et al. 1990 



Bennett et al. 2003al ). Let denote the temperature anisotropy at a sky pixel i. The raw 
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data in a certain band is a set of temperature differences d between pixels in the sky. From 
N observations we have the following observation equations 

t\+ — ti~ = d\ 

t_/V+ — tfq~ = d]\f . 

The above equation system can be expressed by matrix notation 

At = d . (2) 

Where the scan matrix of the experiment A= (a(k,i)), k = 1,---,N and i = 1, ■■■,£/ 
with L being the total number of sky map pixels. The most of elements a(k, 1) = except 
for a(k, i = k + ) = 1 and a(k, i = k~) = —1, where k + denotes the pixel observed by the 
plus-horn and k~ the pixel observed by the minus-horn at an observation k. 

The normal equation of Eq. [T] or Eq. [2] is 

Mt = A T d (3) 
with M = A T A. The least-squares estimate of the sky map results from solving Eq. [3] 

t = M x A T d . 



The WMAP team (IHinshaw et al. 20031 ) uses the following approximate formula to compute 
the iterative solution 

t (n+i) = M- 1 (A T d - A T At (n) ) + t (n) , (4) 

where M -1 = diag(^-, ■ ■ •) is an approximate inverse of M with N{ being the total 
number of observations for pixel i. 

The use of approximate inverse matrix Mr 1 is not necessary. Here we derive an iterative 
formula directly from the normal equation. The Eq. [3] can be expressed as 

fc+=i k~=i k+=i k~=i 

(i = l,2,...,L) . 

Where Ylk+=i means summing over iV^ observations while the pixel i is observed by the 
plus-horn and J2k-=i means summing over N~ observations while the pixel i is observed by 
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the minus-horn, and the total number of observations for the pixel i is Ni — + N i . From 
the above equations we can derive the following iterative formula 

'. ( " +I, = ^(I> +£>)-£ (4-tJ?)), 

fc+=i k =i 

(i = l,2,..-,L) . (5) 

With Eq. HI or Eq. [5] when the number n of iteration is large enough, we get the final 
solution tj = t\ for each pixel i. The Eg. H]used by the WMAP team is an approximate 
formula and Eq. [5] is an exact one, but both has good performance for the differential data 
of a noiseless instrument. With Eq. [5] we can easily study the statistical and systematical 
errors induced by instrument noise, inhomogeneity of sky coverage, inhomogeneity of sky 
temperature, and unbalance between two sky side measurements. 



2. Exposure-dependent noise 



We can directly see from Eq. [5] that their exists exposure dependent noise in a recovered 
temperature map. The WMAP differencing assembly has instrument noise n per observation, 
the real observation data dk = —tp.- +rik. The instrum ent noise is not neg ligible, e.g. the 
mean rms noise o"o = 6.7 mK for the W4-band of WMAP (jLimon et al. 20031 ). The obtained 
temperature differences dk can be taken as a random variable with a standard error cr - The 
final iterative solution ti from Eq. [5] has a noise component being the mean of Ni variables, 
i.e. the temperature ti in a WMAP map for a sky pixel i has an exposure-dependent error 



1 



a 



(6) 



By analyzing CMB maps from the first year WMAP (WMAP1) data, Tegmark et al. 
(2003) find both the CMB quadrupole and octopole having power along a particular spatial 
axis and more wo rks (Ide Oliveira-Costa et al. 2004 ; lEriksen et al. 2004al ; ISchwarz et al. 2004 
Jaffe et al. 20051 ) find that the axis of maximum asymm etry tends to lie c lose to the ecliptic 
axis. A similar anomaly was also found in COBE maps (ICopi et al. 20061 ). The unexplained 
orientation of large-scale patterns of CMB maps in respect to the ecliptic frame is one of 
the biggest surprises in CMB studies (jStarkman and Schwarz 20051 ). A notable asymmetry 
of temperature fluct uation power in two o pposing hemisphere s is also found in the WMAP1 
and COBE results feriksen et al. 2004bl ; Irlansen et al. 2004h . After the release of WMAP 

results in 2006 March, similar large-scale anomalies are still detected in the WMAP3 data 

f lAbramo et al. 20061 : Ijaffe et al. 20061 : ICopi et al. 20071 : lLand & Magueiio 20071 : lEriksen et al. 2007 
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Park et al. 20071 ; IVielva et al. 20071 ; ISamal et al. 20081 ). The unexpected large-scale anoma- 



lies in CMB maps are extensively studied with different techniques, but their reasons still 
remain unclear. Here we show that the exposure-dependent noise should be an important 
source of detected anomalies. 



2.1. Large-scale non-Gaussian modulation 

Because the sky coverage of WMAP mission is inhomogeneous - the number of ob- 
servations bdn^_greates_La^the ecliptic poles and the ecliptic plane being most sparsely 



observed (IHinshaw et al. 20071 ). from Eq. [H]we know that the WMAP temperature maps con- 
tain inhomogeneous exposure-dependent noise components. The remarkable fact, the feature 
of large-scale anomalies detected in WMAP maps being very similar with the WMAP expo- 
sure pattern, strongly indicates the existence of such observational effect on WMAP maps. 
To address the large scale anomalies, such as asymmetry, alignment an d low I power issues 



detected in WMAP data with different techniques, the WMAP team (ISpergel et al. 20061 ) 
describe the observed temperature fluctuations, t, as a Gaussian and isotropic random field, 
t, modulated by a function /(n) 

f(n)=t(n)[l + /(n)] 

where /(n) is an arbitrary modulation function. They expand /(n) in spherical harmonics 



1=1 l=—m 

and use maximum likelihood technique with a Markov Chain Monte Carlo solver to get the 
best fit values of fi m with l max = 2 for the WMAP3 V-band map. The bottom panel of 
Fig. 1 is obtained based on the best fit coefficients, showing in a unifying manner the large 
scale anomalies in WMAP temperature fluctuations which is the same feature that has been 
identified in a number of papers on non-Gaussianity. We calculate the spherical harmonic 
coefficients fi m with l max = 2 for the map oil/N with N being number of observations per 
sky pixel from the WMAP3 V-band data. The top panel of Fig. 1 shows the map of 1/N 
reconstructed based on the coefficients fi m . The middle panel shows the reconstructed result 
for the observation fluctuation map - the map of l/a(N) where the rms variation cr(N) = 
y/((N- (AO) 2 ) calculated within a region of ~ 1° side dimension for each sky pixel. The 
large-scale non-Gaussian modulation features of WMAP temperature map and scan pattern 
being similar for each other suggests that large scale anomalies detected in WMAP maps 
are most probably resulted from observation effect, not cosmological origin. In comparing 




Fig. 1. — Large-scale non-Gaussian modulation features with l max = 2 for the WMAP3 V- 
band data . Top panel: large-scale feature of 1/N map. Middle panel: large-scale feature of 
l/cr(A r ) map. Bottom panel: Best-fit large-scale modulation function f(n) for temperature 
map. 
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Fig. 2. — Latitude distributions from the WMAP3 foreground-cleaned Ql-band map in 
ecliptic coordinates. Black line: rms variation in the map, ctq = 2.27 mK. Green line: 
residual rms after excluding the exposure-dependent anisotropy component by Eq. [71 Red 
line: average number of observations per sky pixel. Blue line: average map temperature (in 
mK, times by a factor of 2.76 x 10 4 ). 



the top and middle panels in Fig. 1, the modulation pattern for the observation fluctuation 
map is more similar to the detected anomalies shown in the bottom panel, indicating that 
the fluctuation of observation numbers could produce additional noise component to the 
recovered temperature map. 



2.2. Alignment and planarity 

To show the anisotropy noise, we compute the average rms variation J (i 2 ) for the 
WMAP3 foreground- cleaned Ql-band temperature map at different latitudes; the results 
are shown by the black line in Fig. 2. The red line in Fig. 2 shows that of average number 
of observations per sky pixel. To modify the exposure-induced rms anisotropy, we use the 
following formula to get the residual variation 

rms = sJ({i-{t)Y)-{ol/N) . (7) 
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Fig. 3. — Latitude distributions of temperature rms variation from WMAP3 foreground- 
cleaned maps in ecliptic coordinates, left: VI band, &q = 3.29 mK. right: Wl band, 
cr = 5.83 mK. 




Fig. 4. — Latitude distributions of temperature rms variation from WMAP foreground- 
cleaned maps in ecliptic coordinates, left: WMAP1 W4 band, o"o = 6.7 mK. right: WMAP3 
W4 band, a = 6.7 mK. 



The residual rms of WMAP3 Ql-band , the green line in Fig. 2, shows that the latitude 
dependences are well excluded by using Eq. [71 or, in other words, the noise per pixel in the 
WMAP map is really exposure-dependent in a way described by Eq. [6j 

Similar features are also found in other waveband maps and the combined frequency 
maps (TOH and ILC maps) as well. For example, Fig. 3 shows the results from WMAP3 
VI and Wl bands. The similarity of structures of latitude distribution of rms variation for 
different bands shown in Fig. 2 and Fig. 3 indicate that they are commonly originated from 
the ununiformity of sky exposure. The average pixel rms is 0.20 mK from N = 1.68 x 10 6 
observations for Wl band, 0.15 mK from 1.12 x 10 6 observations for VI band and 0.13 mK 
from 0.84 x 10 6 observations for Ql ban d. The instrument noise <To for Wl, VI and Ql band 



are 5.85, 3.29 and 2.27 mK respectively ( iLimon et al. 20031 ). We find that the ratios between 
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rmsxyiV of different bands are approximately equal to that between o"o of corresponding 
bands, which is what expected by Eq. [6j Fig. 4 shows that the noise anisotropy for the three 
year WMAP W4-band map is almost the same as that for the first year map, not suppressed 
with data accumulation. 

As demonstrated above, the released WMAP temperature maps contain considerable 
exposure-dependent noise. The remarkable similarity between the large scale non-Gaussian 
modulation feature of WMAP temperature fluctuation and exposure-pattern of WMAP ob- 
servation, as shown in Fig. 1, suggests that large scale anomalies detected in WMAP maps 
are most probably resulted from observation effect, not astrophysical or cosmological origin. 

For foreground removal the WMAP team produce the internal linear combination (ILC) 
map t from the five frequency sky maps tj (i = 1, • • • , 5) by t = ^ tfjtj where the weights 
minimize the variance of final map, Var(t), under the constraint £\ Wi = 1. For the region 
outside of the inner Galactic plane and by a nonlinear search, the weights w are found 
to be 0.109, -0.684, -0.096, 1.921, -0.250 for K, Ka, Q, V, and W bands, respectively 



(iBennett et al. 2003d ). Eriksen et al. (2004b) make minimization of the variance under 



the constraint by means of Lagrange multipliers and obtain the solutions are the inverse 
covariance weights Wi = =^ with Cy being the map-to-map covariance matrix. 

It is needed to inspect the effect of exposure dependent noise in the ILC map, which is 
extensively used in cosmological analysis, although the WMAP team warns against its us e 



for CMB studies because of the complex noise properties of this map (ILimon et al. 20031 ). 
The left panel of Fig. 5 shows the latitude distribution of rms deviation from the WMAP3 
ILC map. There exists a variation component monotonically decreasing from the southern 
pole to the northern pole, expressed roughly by the hand-drawn dotted line in the left panel 
of Fig. 5. The right panel shows the residual rms deviations after subtracting the north-south 
asymmetry component. From Fig. 5 we see that the procedure of linear combination with 
respect to the covariance between different frequencies can effectively depress the instrument 
noise. The latitude dependence of rms variation in the ILC map is weaker, but still visible 
with a structure similar to what observed in W, V and Q band. We can not make modification 
for the ILC map by sky exposure like what we do for Q, V and W data sets, because the ILC 
map is reconstructed by re-adding the data at different frequencies, it is difficult to assign 
an observation number iVj to a pixel i. 
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Fig. 5. — Left: RMS deviation of temperatures vs. ecliptic latitude for the three-year ILC 

map. The vertical coordinate is \J (t 2 ), where t denotes temperature anisotropies (in mK) 
at a certain latitude. Right: Residual rms deviations obtained from the curve shown in the 
left graph by subtracting the north-south asymmetry component (the dotted line in the left 
graph) 



2.3. Pseudo sources 



In studying the nature of discrete 7-ray sources discovered by the COS-B satellite 
experiment, we found that quite a lot apparent discrete sources revealed in the cross- 
correlation map with high significances (about half of sources in 2CG catalog) are pseudo- 
sources produced by th e fluctuation of structured diffuse background of the galactic plane 



( iLi fc Wolfendale 19821 ). Here we use cross-correlation as another indicator to inspect the 
possible large-scale anomaly caused by the unevenly distributed noise in WMAP maps. We 
compute the cross-correlation function c of the W4-band map t and beam profile B by 
c = B T t, and the standard deviation a c of the correlation map. The pixels with |c| > 3cr c 
in ecliptic coordinates are shown in the left panel of Fig. 6. That the apparent features of 
point-like sources concentrated along the ecliptic plane are just generated by observation 
effect as this feature disappears after modifying WMAP temperatures by a factor of 



'(M) 

for both one-year and three-year data, as shown by the right column of Fig. 6. The result 
shown by Fig. 6 remind us that one should take a care before to claim a finding of CMB 
anomaly from observation maps before carefully considering the effect of exposure dependent 
noise. 



Fig. 6. — Cross-correlation map in ecliptic coordinates. Points show pixels with c > 3a c , 
c - cross-correlation of beam profile and sky temperature, a c - standard deviation of c. 
Left panel, from W4-band map. Right panel, from modified W4-band map with multiplying 



A remarkable character in Fig. 5 is the monotonically decline trend of rms variation 
from the south pole to the north pole in the ILC map. To show the north-south asymmetry 
quantitatively, we calculate the asymmetry ratio rms(/3)/rms(— (3) for different latitude f3, 
where rms(/3) is the average rms of sky pixels with latitude within (3 ± 1.5° in the northern 
hemisphere and rms(— j3) is that in the south. The result is shown in Fig. 7. From Fig. 7 
we see that there exists a systematic decrease of rms deviation from south to north with 
maximum asymmetry amplitude ~ 10% — 20%. For different ecliptic latitudes we calculate 
residual rms from the WMAP3 V-band foreground- cleaned map after modifying the exposure 
ununiformity by using Eq. the result is shown in Fig. 8, where we also find a clear north- 
south asymmetry. 

To see if the north-south asymmetry also exists in WM AP exposure, we calculate the 
exposure asymmetry ratio N{(5)/N{—(3) for Wl, VI and Ql band separately, where N{(5) is 
the three-year observation number contributing to the pixels with latitude within (3 ± 1.5° 
in the northern hemisphere and N(—f3) is that in the south, the results are shown in Fig. 9. 
From Fig. 9 we see that there exists north-south asymmetry in WMAP sky survey: exposure 
increases from south to north with maximum amplitude ~ 2% — 4%. It is natural to explain 
the rms asymmetry by the exposure asymmetry as ~ 10% — 20% decreasing of rms deviation 
is just expected by ~ 2% — 4% increasing of exposure, provided that the linear combination 
can not suppress the common and systematic north-south asymmetry in multi frequencies 



temperatures by 




Top panel, from WMAP1 data. Bottom panel, from WMAP3 data. 



2.4. North-south asymmetry 



20 40 60 80 

Latitude /3 (deg) 



Fig. 7. — North-south asymmetry of rms deviation of temperatures in ILC map. (3 - ecliptic 
latitude. rms(/3)/rms(— (3) - ratio between rms deviations at the latitude j3 in the northern 
hemisphere and at that in the southern hemisphere. 



and the exposure dependent noise follows Eq. EJ rms oc 1/ ^/(Ni). 

An apparent asymmetry in the distribution of fluctuation power in two opposing hemi- 
spheres with a remarkable absence of power in the vici nity of the northern ecliptic pole is ob- 



served from the WMAP1 Q, V, and W band sky maps ( jEriksen et al. 2004bl ; lHansen et al. 2004 ) 
By analyzing the issue of power asymmetry with the WMAP3 ILC map and a model of an 
isotropic CMB sky modulated by a dipole field, Eriksen et al. (2007) find that the modu- 
lation amplitude is 11.4% and that the results on hemispherical power asymmetry are not 
sensitive to data set or sky cut. All of these features can be explained naturally by the effect 
of exposure dependent noise as we discuss above in this section. 
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Fig. 8. — Ecliptic latitude dependence of residual temperature rms fluctuations of the 
foreground-cleaned WMAP3 V-band map. The rms fluctuations are modified by using Eq. [7J 
The best-fit line y = -4.46 x 10" 5 a; + 0.085. 




Fig. 9. — North-south asymmetry of WMAP3 exposure. N(/3) is the total number of obser- 
vations contributing to the pixels with ecliptic latitude from (/3 — 1.5) deg to (/3 + 1.5) deg in 
the northern hemisphere, N(—(3) is that in the southern hemisphere. Left panel, Wl band. 
Middle panel, VI band. Right panel, Ql band. 
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3. Foreground contamination 

Besides the statistical noise, another kind of error in recovered WMAP CMB maps is 
systematic distortion. From Eq. [5] we can see that in an iteration the temperature estimation 
at a sky pixel i is evaluated based on temperatures of many pixels on a circle in the sky 
sphere beam (beam separation) apart from the pixel i, including N~ pixels pointed by minus- 
horn when plus-horn pointing to i and pixels pointed by plus-horn when minus-horn 
pointing to i. The sky temperature tk+ observed by plus- horn and by minus-horn are 
differently placed in the right side of the iterative formula Eq. [51 The iterative solution t« 
could be deviated from the true temperature ti due to inhomogeneity of temperature sky 
and unbalance between two sky side beam measurements. 



3.1. Map distortion by a hot source 

The foreground induced systematic effect on a recovered map is not only limited in the 
region containing foreground sources, but spreated over the whole sky. A hot foreground 
source pointed by the side beam A of a radiometer can distort the recovered temperatures 
of sky pixels pointed by the side beam B with separation angle 9b eam to the source through 
map-making iterations. From Eq. [5] (or Eq. HJ) the first iterative solution of temperature of 
a sky pixel i can be expressed as 

t{ ti (t t^ ^ring (8) 

where ( ) r i ng denotes averaging on the scan ring with separation angle 6b eam to the pixel i. 
For zero initials t (0) = 

ti (fyring i (9) 

a hot source contained on the scan ring will let (t) r i ng > and the recovered temperature 
t!p < 0, indicating that a hot foreground source might systematically make the recovered 
temperatures on its scan ring lower. 

A one-dimensional simulation is made to show such effect. A true temperature "sky" 
containing 500 pixels is produced as a white noise series of zero mean and 0.2 mK standard 
deviation with a hot source in the pixel interval 240-260, as shown in the upper panel of 
Fig. 10. To produce differential data, we simulate a scan process of a differential radiometer 
with instrument noise <jq = 4.0 mK and beam separation 6b ea m — 100 pixel. From i = 1 
to 500, for each pixel % we produce 2500 temperature differences £j — tj + n where j = 
i + 100 for i < 400 or otherwise j — i — 400, n is sampled from the normal distribution 
with zero mean and standard deviation ctq. We use Eq. H]to reconstruct temperature map 
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from the differential data. For pixel k = 350 the recovered temperatures from first fifteen 
iterations, t k , t£\ • • • , t^ 5 " 1 are shown in Fig. 11, where we see that the recovered temperature 
of first iteration for the pixel 9f, earn away from the hot pixel really drops down dramatically 
as expected by Eq. M and after a few iterations converged to a stable value which still 
significantly lower than its real temperature. The lower panel of Fig. 10 shows the recovered 
temperature map after 50 iterations. There exist two significant cold regions in the recovered 
map with pixels 340 — 360 and 140 — 160, each have a distance of 6^ eam to the hot source 
interval and their temperature variation negatively related to that of the hot source. Two 
more regions distorted in the similar way can be found in the recovered map, each placed 
100 pixel away from the above two cold regions respectively, pixel intervals 440 — 460 and 
40 — 60. The above simulation result indicates that a hot source can distort temperatures 
in the map recovered from differential data at pixels away from the source with distance of 
the beam separation and such distortion can propagate to more far away pixels. 
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Fig. 10. — True and recovered temperatures. Upper panel: True temperatures. Lower panel: 
Recovered temperatures. 
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Fig. 11. — Recovered temperature at pixel 350 vs. step number of iteration by map-making 
Eq. HI The one-dimensional true temperature map, shown in the upper panel of Fig. 10, has 
a hot source with a peak at pixel 250, and the differential data is simulated for a radiometer 
with beam separation of 100 pixel. 



3.2. 141° rings in WMAP maps 

The beam separation angle 9 beam of WMAP radiometers is 141°. We predict from the 
above simulation that there should exist strongest negative correlation in WMAP tempera- 
ture maps between temperatures of two sky pixels separated 141° each other. This prediction 
is confirmed by analyzing temperature maps in Q, V and W bands released by WMAP team. 
For the WMAP3 Q -, V- and W-band maps with HEALPix resolution parameter iV^e — 128 



( jGorski et al. 20051 ). we calculate cross-correlation coefficients C ta (T b ){9) between t a and (t b ) 
for different separation angle 9, where t a denotes the temperature of a sky pixel a and (tb) 
the average temperature of the ring with separation angle 9 to the pixel a, the obtained 
distributions are shown in Fig. 12. 

From Fig. 12 we find that the strongest negative correlation appear indeed around 141° 
separation in the correlation distributi on for each band. If p ixel a is limited within the sky 



region of the foreground mask Kpl2 (IBennett et al. 2003bl ) which consists mainly of hot 



pixels, around 141° separation show more strong negative correlation just as expected from 
assuming that hot foreground emission induce the correlation. 

The program synfast in HEALPix software package (available at jhttp:/ /healpix.jpl.nasa.gov 
can create temperature maps computed as realizations of random Gaussian fields on a sphere 
characterized by the user provided spherical harmonic coefficients of a angular power spec- 
trum. To test the significance of the correlation C(141°) = —0.238 in Q-band WMAP map, 
simulated temperature maps are created with the synfast program from the best fit ACDM 
model power spectrum with Q-band beam function and noise property. For each simulated 
CMB map, we compute the correlation coefficient between temperatures of the pixels in 
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Fig. 12. — Correlation coefficient C ta (t b )(Q) between temperatures t a and (tb) vs. separation 
angle 6. t a denotes the temperature of sky pixel a, (tb) the average temperature of the ring 
with separation angle 9 to a. Filled square: pixel a only within the foreground mask Kpl2. 
Square: for whole sky. Triangle: only for the sky region out of the mask KpO. Top panel: 
WMAP3 Q-band. Middle panel: WMAP3 V-band. Bottom panel: WMAP3 W-band. 



Kpl2 region and average temperatures on their 141° rings. From 1000 simulated maps we 
get (7(141°) = 0.108±0.080, indicating that the 141° negative correlation detected in WMAP 
Q-band map has a significance of 4.3cr. Similarly, for V-band we get C(141°) = —0.231 with 
significance 4.4<r evaluated by C"(141°) = 0.099 ± 0.075 from simulated CMB maps, and for 
W-band C(141°) = -0.268 with significance 5.0a from C"(141°) = 0.094 ± 0.072. 

It is expected from our simulation shown in Fig. 10 that foreground hot sources can 
produce cold rings of 141° separation in WMAP temperatures. To test this effect, for 
each sky pixel we calculate the average temperature of its 141° ring from the WMAP3 
Q-band, V-band, W- band and ILC temperature maps with N si d e = 128 and KpO mask 



( iBennett et al. 2003bl ) for foreground clean. The result is shown in Fig. 13. There is a low 




Fig. 13. — Average temperature maps of 141° ring in Galactic coordinates from WMAP3 
Q-band, V-band, W-band and ILC map (from top to bottom panel). 
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temperature region near Galaxy center in Fig. 13 for each band indicating that most of the 
141° rings corresponding to the galaxy center region are cold. This is consistent with our 
expectation. 

To roughly estimate the magnitude of temperature distortion by foreground emission, 
we pick up 2000 hottest pixels in the 3-year WMAP Q-band map and then find out their 141° 
rings. The average temperature of all pixels on these rings and out of the KpO mask (cover 
about 14% of the sky, see Fig. 14) is calculated to be (t) = —0.0117 mK. The correspondent 
value from 1000 simulated CMB maps is (t 1 ) = 3.4 x 10~ 5 ±0.006 mK. Similarly, for V-band we 
get (t) = -0.0125 mK and (if) = -4.7 x 10~ 5 ±0.006 mK, for W-band (t) = -0.0141 mK and 
(£') = 8.8 x 10 -5 ± 0.006 mK. Therefore, in foreground-cleaned WMAP maps, considerable 
foreground-induced deviations still exist with amplitude comparable to the fluctuation of 
CMB signal and over broad sky region. 




Fig. 14. — 141° rings of 2000 hottest pixels in 3-year Q-band WMAP map, after using the 
KpO mask, in Galactic coordinates. 

The foreground- induced deviation in WMAP temperatures should also distort the tem- 
perature angular power spectrum. We use simulated temperature maps created with the 
synfast program to roughly estimate the magnitude of foreground-induced distortion for the 
WMAP power spectrum on large-scale (low-/) region. For each simulated map, we sub- 
tract 0.01 mK from each temperature of pixels on the 141° rings corresponding to the 2000 
hottest pixels (shown in Fig. 14) to get a distorted map. We compute the power spectrum 
for simulated and distorted map respectively and find that, on average, the deviation be- 
tween the original and distorted spectrum densities at / = 2, 3 is ~ 10%, indicating that the 
foreground-induced distortion on WMAP power spectrum cannot be ignored for a precision 
cosmology study. Strongest hot sources may produce cold spots by the foreground induced 
distortion out of the foreground region. Temperature distortions in a CMB map caused by 
foreground sources with different scales and observation dependent noise could distort the 



-19- 



angular power spectrum on wide range of angular scale. 



3.3. Large non-Gaussian spots 



A large cold spots with a radius of ~ 5° centered at (I, b) = (209°, —57°) and the 
lowest temperature ~ —0.15 mK h as been detected i n WMAP 1 and WMAP3 maps with 



the wavelet and other techniques (Vielva et al. 20041 ; ICruz et al. 20051 ; ICruz et al. 2006 



Cruz et al. 2007al ; IVielva et al. 20071 ). and was explained as a cosmic texture, a remnant of 
symmetry-breaking at energies close to the Planck scale in the very early universe (ICruz et al. 2007bl ) 
We find that the 141° rings of the pixels in the spot region across through the hot Galactic 
plane with temperature up to ~ 5 mK, therefore, the foreground induced systematic effect on 
recovered maps discussed in §3.1 should be a more plausible explanation to the detected large 
cold spot than the texture hypothesis. The foreground contamination and other systematic 
effect may also produce other more detected large spots in WMAP maps. 



4. Discussion 

4.1. Can the exposure induced anisotropy be corrected ? 

It has to be pointed out that Eq. [7]can be used to modify the effect of exposure dependent 
noise, like what we do for the latitude distribution of rms variation shown in Fig. 2, only for 
the case that the map rms fluctuation (t 2 ) itself is the directly analyzed quantity. However, 
almost all analysis works are based on CMB maps of temperature. We know from Figs. 2- 
5 that the released WMAP temperature maps contain considerable exposure dependent 
noise. It is no way from a recovered temperature map to produce a corrected map in which 
the instrument- induced and exposure-dependent noise can be eliminated. Therefore, the 
anisotropy noise should contribute to large scale anomalies detected in existent CMB maps. 

The exposure inhomogeneity of WMAP comes from its scan strategy, which can't be 
suppressed through accumulating observation time. To avoid the observation effect, it is 
needed to remake temperature maps from a uniform differential data set obtained by giving 
up partial observation data for pixels of high exposure. Comparing released WMAP maps 
and new maps from uniform data will help us to judge their origin of detected large scale 
anomalies, e.g. the low I power issues detected in WMAP data, the unexplained orientation of 
large-scale patterns of CMB maps in respect to the ecliptic frame, the north-south asymmetry 
of temperature fluctuation power etc., and to see if the observational effect can also influence 
the angular power spectrum as well. 
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4.2. Avoiding the foreground-induced distortion 



The distortion by hot foreground sources on their 141° rings in a WMAP map can not 
be removed with a foreground mask on the recovered temperature map. What's needed is 
to use the mask on the original differential data before map-making to avoid the foreground- 
induced error in the recovered map. The top panel of Fig. 15 shows the difference between the 
recovered and true temperature distributions (shown in the lower and upper panel of Fig. 10 
respectively), where the distortion structure caused by the hot source on pixel 240 - 260 and 
beam separation of 100 pixel is clearly shown. We redo the temperature reconstruction with 
excluding the temperature differences that contain a beam side pointing to a pixel between 
240 - 260 ("mask region") during iterations for the pixels out of mask, the result is shown 
in the middle panel of Fig. 14, where the distortion structures are really suppressed. 

A weakness of using mask in map-making process is decreasing the number of useful 
differential data. Another approach to avoid the distortion in recovered map by foreground 
emission is to properly set the initials of iteration for the foreground region. From Eq. [S] we 
see that the temperature deviation of the first iterative solution, t\ — U, will be suppressed 
if the temperature initials at pixels of hot source are set to be close to their true values to 
let (t - t (0) )rin 9 « 0. The initial tf ] of pixel % can be taken as 



Where Ylk+=i means summing over the observations while the pixel i is observed by the 
plus-horn and the pixel pointed by the minus- horn is out of mask, ^ fc _ =i means summing 
over the observations while the pixel % is observed by the minus-horn and the pixel pointed by 
the plus-horn is out of mask, N[ is the total number of used observations. For the simulated 
differential data from the true temperatures shown in the upper panel of Fig. 10, we make 50 
iterations with Eq. H] starting from initials calculated by Eq. [TUJ the distortion structures are 
satisfactory suppressed in the resultant solution, as shown in the bottom panel of Fig. 15. 



We demonstrate in this paper that for existent CMB maps both the observation depen- 
dent noise and systematic error induced by foreground emission can not be neglected and 
both can produce large-scale anomalies and distort the angular power spectrum. These errors 
can not be completely excluded by performing noise suppressing or using foreground mask 
on temperature maps. We suggest to remake temperature maps from the original WMAP 
time-order-data by a modified algorithm with applying foreground mask in map-making to 




(10) 



4.3. Remaking WMAP maps 
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Fig. 15. — Differences between recovered t and true temperature t . The one-dimensional 
temperature distribution to consists of white noise and hot source between 240 - 260 pixel 
(mask region), as shown in the upper panel of Fig. 10. Differential data is obtained by 
simulation of one-dimensional scan for to by a differential radiometer with beam-separation of 
100 pixel. Recovered t are calculated with map-making Eq. HJafter 50 iterations. Upper panel: 
t recovered from all differential data (the lower panel of Fig. 10). Middle panel: t recovered 
by iterations from all data for mask region and with excluding differences that contain hot 
source temperature for pixels out of mask. Bottom panel: t recovered by iterations with 
initials estimated from the differential data by Eq. [TUJ 
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exclude mask pixels from use in iterations for CMB dominated region (or properly set tem- 
perature initials before iteration), and/or keeping used differential data uniform by giving up 
partial observation data for pixels of high exposure. New maps from modified map-making 
algorithm will help us to judge the origin of large scale anomalies detected in released WMAP 
maps, e.g. the low / power issues, the unexplained orientation of large-scale patterns in re- 
spect to the ecliptic frame, the north-south asymmetry of temperature fluctuation power and 
the large non-Gaussian spots, and to see to what extent the statistical and systematical errors 
influence the angular power spectrum and the derived cosmological parameters. Believable 
conclusions on CMB anisotropy anomalies and precise temperature angular power spectrum 
from differential measurement should be based on temperature maps with homogeneous sky 
exposure and should avoid foreground-induced distortion during map-making. 

This study is supported by the National Natural Science Foundation of China and the 
CAS project KJCX2-YW-T03. The data analysis in this work made use of the WMAP data 
archive and the HEALPIX software package. 
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